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Summary 

This paper analytically examines the unsteady fluid dynamics of a 
vortex filament subsequent to a normal collision of the vortex with a 
solid body. In particular, the breakdown or reconnection 
phenomena, post-collision, for a vortex filament is studied. The 
paper does not investigate the collision dynamics process itself. The 
derived exact solution is based upon the laminar viscous form of the 
Helmholtz equations. 


Nomenclature 

r Radial coordinate, origin at filament axis, ft 
r* Nondimensional radial coordinate, r*=r/s 
r c Vortex filament core size radii as a function of time, t, ft 
Re Vortex filament Reynolds number, Re = y/v 
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Axial distance of filament breakpoint (time equal zero) from 
origin, ft 

Time, sec 

Nondimensional time parameter, t*=vt/s 2 

Velocity vector, cylindrical coordinates, V = [v r v e vj, ft/sec 

Radial velocity component, ft/sec 

Nondimensional radial velocity component 

Axial (along vortex filament axis) velocity component, ft/sec 

Nondimensional axial velocity component 

Tangential velocity component, ft/sec 

Nondimensional tangential velocity component 

Axial (along vortex filament axis) coordinate, origin at 
intersection of filament segments' plane of symmetry and 
filament axis, ft 

Nondimensional axial coordinate, z*=z/s 

Vortex filament initial circulation strength, ft 2 /sec 

Vorex circulation, ft 2 /sec 

Kinematic viscosity, ft 2 /sec 

Angular coordinate, radians 

Vorticity vector, co = [co r to e ©J , ft 2 /sec 


Quasi-stream function 


Introduction 


Vortex filament 'collisions' with solid bodies occur in a number of 
real world examples for rotary-wing aircraft: blade vortex 
interactions where vortices shed from one blade collide with another 
blade; trailed tip vortex interactions/collisions with a helicopter 
airframe/tailboom in hover and low-speed flight; or, alternatively, 
proprotor trailed tip vortex collisions with wing surfaces for tiltrotor 
aircraft in airplane-mode. For, blade vortex interactions, references 
1-3 provide qualitative discussion regarding the vortex filament 
collision process. References 4-7 provide a similar qualitative 
discussion of vortex fuselage/airframe/wing collisions. 

Very little work to date has been performed examining the post- 
collision dynamics of vortex/solid body interactions. Most work has 
either concentrated on the collision process itself (see reference 8) or, 
alternatively, has looked at vortex-on-vortex collisions (see reference 
9, for example). A related academic exercise is the study of ring 
vortices throughout the collision process; see reference 10, for 
example. 

Nonetheless, flow visualization techniques are sufficiently maturing 
so as to examine in detail the rotary-wing/airframe vortex filament 
collision process - see, in particular, references 4, 5, and 7. It's clear 
that current rotor wake CFD analyses are inadequate to represent 
rotor vortex filament collision and subsequent breakdown or 
reconnection. This paper proposes as a initial step towards 
understanding the very complex flow phenomena of rotor wakes a 
first order analytical treatment of the vortex filament/solid body 
collision process. This paper will concentrate on the post-collision 
vortex filament fluid dynamics. 
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Problem Description 


Figure 1 illustrates the point that there is going to be varying 
degrees of vortex filament breakage (in form of filament distortion 
and separation of breakpoints) depending on the severity of the 
collision process. Vortex filament collisions with thin solid bodies 
will tend to minimize the separation distance between the filament 
segment breakpoints as well as minimize the filament distortion. 
Therefore, collision with thin solid bodies will tend to result in the 
'reconnection' of the vortex filaments. Conversely, collisions with 
bluff bodies will tend to result in large breakpoint separation 
distances and considerable filament segment distortion - therefore, 
most likely resulting in filament segment 'breakdown.' 


V 






Collision with Thin Body Collision with Bluff Body 


Figure 1 - Hypothetical Vortex Filament Distortion/Breakage After 

Collision with Solid Bodies 
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Figure 2 illustrates an idealized version of the vortex filament in the 
intermediate process of vortex reconnection/breakdown, post- 
collision. Distortion of the filament as to deflection of filament axis 
laterally from its undisturbed orientation is neglected. Initial axial 
flow (along the filament) is assumed to be zero as a result of the 
collision process: the filament, as it wraps around the solid body, will 
have to obey the no-flow boundary constraint. 



Vortex 

Filament 

Segment 


Fi gure 2 - Vortex Geometry Description. tVortex Filament Described 
by Contour of Maximum Velocity") 
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Governing Equations 


The Helmholtz equation for laminar-flow vorticity (references 13 and 
14) is given by 

— — = (to • V)V + vV 2 co (1) 


The term Dco/Dt is the 'particle derivative' for the vorticity. And, 
further, (references 13 and 14) 


divco = V • co = 0 


( 2 ) 


where 


co = curlV = V x V 


(3) 


Considering the vorticity component about the z-axis (cd 2 ), for the 
problem under discussion — and using the cylindrical coordinate 
system for the Helmholtz, vorticity divergence equations, and the 
vorticity definition -- the above two equations become respectively 
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(5) 
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Imposing 9-symmetry on equations 4, 5, and 6a-c gives 
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(9a-c) 


0 ) 
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Note that substitution of equation 9a into equation 8 automatically 
satisfies equation 9c. Finally, the continuity equation (0-symmetry) 
requires that 


10/ v 1 0V„ 0V. v r 0V r 0V, 

(rvJ + --^ + -^^^ + ^ + —L = 0 
r or r 30 3z r dr dz 


( 10 ) 


Or, rather 



This completes the set of required governing equations necessary to 
solve the described problem for laminar flow. 


Solution 


A class of unsteady laminar flow problems called 'moving boundary' 
problems has been exhaustively studied in the literature. By 
assuming that v r =v z =0, the Navier-Stokes equations (or, 
correspondingly, the Helmholtz vorticity equation) can be reduced to 
the well-known unsteady heat conduction equation (reference 15). 
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The Lamb-Oseen vortex solution (references 13 and 14) is one 
member of this class of unsteady laminar flow problems. 

What is sought then is a new class of flow problems that not only 
encompasses the 'moving boundary' class of problems but the 
problem described in this paper. But key to defining this new class 
of flow problems remains the mathematical challenge of reducing the 
Helmholtz equation, with 0-symmetry, to the unsteady heat 
conduction equation. 

It is proposed that this new class of solution can be defined by 
assuming that the radial velocity component, v r , is not merely equal 
to zero but instead is proportional to the tangential velocity gradient 
with respect to the z-axis (equation 12). (Note that if the tangential 
velocity gradient is zero, such as in the case of the Lamb-Oseen 
vortex, then the 'moving boundary' class of problems is recaptured.) 


fK 

dz 


Or 


.3V,, 


(12a) 


(12b) 


Where the length scale factor, l, is a constant that transforms the 
proportional relationship of equation 12a to the equivalence 
relationship of 12b. (The constant l has the unit of length - hence it 
being called a 'length scale factor.') This will be discussed in detail 
later in the paper. 

Meanwhile, continuing with the analysis with the proportional 
relationship for radial velocity, applying equation 12a to equation 11 
gives a second proportional relationship based on flow continuity for 
the axial velocity. 
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(13) 
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Applying equations 12a and 13 to equations 9a-c gives 


<D r - _v r 
<°z °c -v 2 


( 1 4a-b) 


Application of equations 12a, 13, and 14a-b to equation 7 causes the 
convective acceleration terms to cancel out the vortex stretching 
terms and thereby achieves the objective of reducing the Helmholtz 
equation to the unsteady heat conduction equation (equation 15). 


3t °z . Igggt . 

dt \ dr 2 r 3r dz 2 \ 


(15) 


Reference 15 reveals that a particular solution of equation 15, for a 
instantaneous point source, is given by 


to. 


point source 


C -(r 2 +(z-z.) 2 )/4vt 

8(jtvt) w 5 


(16) 


where c is a constant, z* is the origin of the instantaneous point 
source with respect to the z-axis, and r and z are the coordinates of 
the point for which the vorticity is predicted. 

Correspondingly, reference 15 derives the solution to the unsteady 
heat conduction equation for an instantaneous line source of infinite 
length (represented as a continuous distribution of instantaneous 
point sources) by the equation 17. 
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CO 


z linfinite line source 
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47CVt 


( 17 ) 


With the constant, c, set equal to the circulation strength, y, 

(reference 14), then equation 17 becomes the vorticity distribution 
for a Lamb-Oseen vortex. 

To solve the unsteady laminar flow problem for a rectilinear vortex 
filament that has been 'cut' or ’broken,' however, it is necessary to 
represent the vortex filament segments with semi-infinite 
distributions of instantaneous point sources (equation 18). Note that 
the constant, c, still has been assigned the value c=y. 
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(18) 


The integral terms in equation 18 can be found in standard 
handbooks for integration formulas (see reference 16). The solution 
of vorticity component with respect to the z-axis is given by 


co z =-^-e' r2/4vt < 

87Wt 


erfc 
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V4vt 
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+ erfc 


_±£.YI 

V4vt )\ 


(19) 


where erfc(x) is the complementary error function, i.e.. 


er fc(x) = -|j ) 


du 


( 20 ) 


Noting that complementary error function is related to the error 
function by the relationship erfc(x)= l-erf(x). And, further, 
erf(-x) = -erf (x) , erf(°°) = i, and erf(0) = 0. With these properties of the 
error function and complementary error function in mind, it can be 



readily seen that equation 19 satisfies three important boundary 
constraints and/or conditions. 


0). -» 


_Y_ e -r 2 /4vt 
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for s — » 0 


Z — > ±oo 


t — > 0 for |z|>s 

(Vorticity distribution approaches that of the Lamb-Oseen Vortex) 


Having derived the vorticity distribution, the axial velocity 
distribution follows straghtforwardly from equation 14b, given 
equation 19. 


— c~ r */ 4vl 

87tvt 




+ erfc 



(21) 


Haivng established the proportional relationship v z °c f(r,z,t), it is 
necessary to multiply equation 21 by a length scale to arrive at an 
equivalence relationship v z =/f(r,z,t). The appropriate length scale 
factor is the separation distance, i.e., ^ = ±s. Further, specifically, 

^ = +s when z<0 and / = — s when z>0. Besides being a natural length 
scale factor, using s as the length scale automatically satisifies the 
'moving boundary' extreme of the vortex filament 
breakdown/reconnection problem; i.e., when s = 0 then v z =0. 


v. = 


_T^_g-r 2 /4vt 

87tvt 



f 

+ erfc 

\ 



where z < 0 


(22a) 
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(22b) 



and where z>0. 


Given equations 9c and 19, the tangential velocity distribution can be 
derived through solution of a first-order ordinary differential 
equation. 
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e — 



-e 



+ erfc 



(23) 


The analogous attributes of equation 23 with respect to the Lamb- 
Oseen vortex are obvious (reference 14). As demonstrated in a 
similar manner with respect to the vorticity distribution, if s -» 0 
then, correspondingly, the tangential velocity distribution reduces to 
the Lamb-Oseen profile. 

Next, given equations 23 and 12a, the radial velocity, v r , can be 
derived by taking the derivative of v e with respect to z (equation 
24). 


2v/rr 


vt • r 


(l-e' rJ/4vr ){ 
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(24) 


Applying the length scale factor, f = ±s (where ^ = +s when z<0 and 
i — — s when z>0), to equation 24, given equation 12a, gives the 
following expression for the radial velocity distribution 


v r =+ = (i- e - f2/,vl ){e 

2Vtt 3 vt • r ' 
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(25a) 


when z<0. 
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(25b) 



when z > 0. A key outcome of equations 25a-b is the prediction that 
the radial velocity is negative (inward to the vortex core) for all 
values of r, z, and t. 

The process and rate by which the vortex filament segments 
'reconnect' can be quantified by a single parameter, % • The 
derivation of this parameter will now proceed. First, the vortex 
segment circulation will be defined by equation 26. 



0 


Substitution of equation 19 into 26 and performing the required 
integration yields. 



(27) 


See reference 14 for the derivation for the circulation of an infinite 
line vortex. What is needed to examine the rate of filament segment 
reconnection is the rate of increase of the circulation in the 
intermediate region (z=0) between filament segments versus the 
circulation of the filament far away from the breakpoints (z — ) ±°°) . 
Therefore, the parameter, %, will be defined as follows 


1 = 


EU 

rU. 


r— >oo 


(28) 


Application of equation 27 to equation 28 — taking the appropriate 
limits for each term in the numerator and denominator of equation 
28 — gives the following expression for %. 



s 


(29) 


X = erfc 


V4vt 


The implications of equation 29 and equations 22a-b, 23, and 25 of 
the unsteady flow behavior of the vortex sements will be discussed 
in the next section of the paper. 


Discussion of Results 


A set of exact solutions have been derived for the unsteady laminar 
flow behavior of a vortex filament subsequent to a perpendicular 
collision with a solid body. A number of simplifying assumptions 
was required in order to derive this solution. Chief among those 
assumptions are: (1) that the filament segments still remained 

straight, parallel, and aligned with each other; (2) a clean 'break' or 
'cut' existed such that no initial vortical flow existed in the 
intermediate region between the breakpoints; (3) implied, but not 
explicitly assumed, was that the vortex filaments could be 
represented by a Lamb-Oseen vortex as an initial condition. 

Given this set of exact solutions, it is necessary to examine their 
implications for the vortex filament reconnection (or breakdown) 
process. In order to examine the results of the above analytical 
treatment the following dimensionless parameters will need to be 
defined: r‘=r/s; z*=z/s; t*=vt/s 2 ; v*=v r (s/y); v*=v z (s/y); v*=v e (s/y). 

Figure 3 is the tangential velocity profile with respect tor r* for 
various values of t*. The time dependent behavior manifested in 
figure 3 not only aacounts for the Lamb-Oseen vortex viscous 
diffusion effects but also accounts for viscous vorticity transport 
from the vortex filament segments (|z|>s) to the intermediate region 
(|z|<s) between the filament breakpoints. 
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* 

r 


Figure 3 - Tangentia l Velocity Profile as a Function of Time at the 
Vortex Filament Breakpoint (z*=-l) 


The normalized tangential velocity profile with respect to r* is similar 
to the well-recognized Lamb-Oseen tangential vortex profile - refer 
to figure 4 - except that the viscous vorticity transport mechanism 
results in a net reduction in magnitude of the tangential velocities. 
There are more significant differences in the flow behavior of the 
reconnecting vortex filament segment problem ~ as will be seen. 
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Figure 4 - Differences between the Lamb-Oseen Vortex Tangential 
Velocity Profile and the Profile for Broken but 'Reconnecting' 
Filaments fz*=-1.0 and t * =1.0) 


As noted earlier, with the Lamb-Oseen vortex , v r =v 2 =0. The radial 
and axial velocities for the 'broken' vortex filament segments are, on 
the other hand, substantial. Figures 5 and 6a-b illustrate the vortex 
filament unsteady flow behavior. Figure 5 shows the vector flow 
field for the r-z plane. 
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Figure 6a-b are streakline plots for the r-0 and r-z plane. The key 
feature of the 'broken' vortex filament unsteady flow are spiral 
streaklines. 





The dependence of the parameter x (equation 29) -- the ratio of 
circulation of the flow at z=0 versus that at z-»±« — with respect to 
the nondimensionalized time parameter, t* is shown in figure 7. 





Interpretation of figure 7 suggests that the greater than vortex 
filament breakpoint separation, s, or the lower the kinematic 
visocosity, v, the slower (or less likely) the filament reconnection will 
be accomplished. 



* 

t 

Figure 7 - v as a function of t * 


There can be no absolute means by which to define when a vortex 
filament will 'breakdown' versus 'reconnect.' If one were willing to 
wait an infinite amount of time all broken or cut vortex filaments 
will eventually reconnect, given one interpretation of equation 29. 
However, as a practical matter, there will quickly reach a point 
where the vortex filament segments will have decayed sufficiently 
via viscous diffusion such that no discernable flow structure will 
exist by the time reconnection is achieved — therefore an observer 
would likely conclude that the vortex filaments have suffered 
breakdown rather than reconnection. Only some sort of semi- 
arbitrary criteria can be established to distinquish between vortex 
filament segment reconnection or breakdown. The author will not 
propose such a criteria in this paper. 


Conclusions 




An exact solution for the breakdown or reconnection of a laminar 
vortex filament — having been 'cut* or 'broken' subsequent to a 
collision with a solid body — has been developed. The analytical 
solution is based upon the definition of a general class of flow 
problems that encompasses as a subset the extensively studied 
'moving boundary' laminar flow problems. The viscous Helmholtz 
equation for vorticity is reduced to the unsteady heat conduction 
equation. Consequently, instantaneous lines sources for the vorticity 
can be used to model the vortex filament breakdown or reconnection 
process. The analytical solution does not treat, though, the actual 
collision process itself. 
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